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ABSTRACT 



In a recent paper (|Moskalenko et all |2002), it has been shown that the flux of secondary cosmic ray (CR) 
antiprotons appears to be contradictory to measurements of secondary to primary nuclei ratios in cosmic 
rays when calculated in the same Galactic propagation model. The contradiction appears as a value of the 
diffusion coefficient necessary to match the secondary ratios p/p and B/C. In particlular, it was shown that 
the reacceleration models designed to match secondary to primary nuclei ratios produce too few antiprotons. 
It is, however, clear that some reacceleration is unavoidable in the turbulent interstellar medium. Here we 
discuss an idea of how to improve reacceleration model by allowing for the damping of interstellar turbulences 
on the small scale by cosmic rays, mostly protons. This would lead to increase in the mean free path lenghts 
at low energies, the well-known phenomena empirically discovered in the Leaky Box models, thus producing 
less secondary nuclei. Antiprotons will remain almost non-affected due to their high energy threshold of 
production cross section. 



INTRODUCTION 

The spectrum and origin o f antiprotons in CR has been a matter of act ive debate since the first reported 
detections in balloon nights (|Golden et all Il97fll : IBogomolov et all Il979h . There is a consensus that most 
of the CR antiprotons observed n ear the Earth are "se condaries" produced in collisions of energetic CR 



particles with interstellar gas (e.g.. iMitchell et al 



The spectrum of secondary antiprotons has a peak at about 2 GeV decreasing sharply towards lower ener- 
gies. This unique shape distinguishes antiprotons fro m other CR species . Over the las t few years the accu 
racy has been improved sufficiently (BESS 1995-2000, Orito et al. . 2000l : Sanuki et al. . 2000l : Asaoka et al 



2002) that we can restrict the spectrum of the secondary component accurately enough to test Galactic CR 
propagation models, and the h eliospheric modulation. 

It has been recently shown (jMoskalenko et ail l200ll . I2OO2I : see also Molnar and Simonl l200l|_ Sina et al 



2002 ) that accurate antiproton measurements during the last solar minimum 1995-1997 (BESS. lOrito et all 
2000) are inconsistent with existing propagation models at the ~ 40% level at about 2 GeV while the stated 
measurement uncertainties in this energy range are now ~ 20%. The conventional models based on local 



CR measurements, simple energy dependence of the diffusion coefficient, and uniform CR source spectra 
throughout the Galaxy fail to reproduce simultaneously both the secondary to primary nuclei ratio and 
antiproton flux. 

The reacceleration model designed to match secondary to primary nuclei ratios, e.g., boron/carbon 
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( Simon et all Il98fit ISeo and Ptuskinl Il994h aroduce too few antiprotons because, e.g., matching the B/C 
ratio at all energies requires the diffusion coefficient to be too large. The models without reacceleration can 
reproduce the antiproton flux, however they fall short of explaining the low-energy decrease in the secondary 
to primary nuclei ratio. To be consistent with both, the introduction of breaks in the diffusion coefficient 
and the injection spectrum is required, which would suggest new phenomena in particle acceleration and 
propagation. 

This forces us to developing a more sophisticated treatment of wave-particle interaction in a course 
of CR propagation in the turbulent interstellar medium. (An alternativ e idea of a local "unprocessed" 
nuclei component in low-energy CR is evaluated in Moskalenko et all 2003.) In particular, we suggest a new 
approach which includes the interaction of particles with waves in the interstellar medium in a self-consistent 
way and take into account the wave damping on energetic particles. This requires numerical methods and 
iterative procedure to derive the diffusion coefficient at arbitrary energy. Secondary-to-primary elements 
ratios and antiproton flux should be used as final indicators of the consistency. In our calculations we use 
CR propagation code GALPROP. This work is currently in progress. 

BASIC FEATURES OF THE GALPROP MODELS 

Th e cylindrically symmetric GALPROP models have been described in detail elsewhere ( Strong and Moskalenko) . 
1998); here we summarize their basic features. 

The models are three dimensional with cylindrical symmetry in the Galaxy, and the basic coordinates 
are (R,z,p) where R is Galactocentric radius, z is the distance from the Galactic plane and p is the total 
particle momentum. In the models the propagation region is bounded by R = Rh, z = ±Zh beyond which 
free escape is assumed. 

The propagation equation we use for all CR species is written in the form: 
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where ip = ip(r,p,t) is the density per unit of total particle momentum, ip(p)dp = 47rp 2 /(p) in terms 
of phase-space density /(p), q(r,p) is the source term, D xx is the spatial diffusion coefficient, V is the 
convection velocity, reacceleration is described as diffusion in momentum space and is determined by the 
coefficient D pp , p = dp/dt is the momentum loss rate, Tf is the time scale for fragmentation, and r r 
is the time scale for radioact i ve de cay. The numerical solution of the transport equation is based on a 
Crank-Nicholson ([Press et al 1 11 9921 ) implicit second-order scheme. The three spatial boundary conditions 
ip(Rh, z,p) = tp(R, ±Zh,p) = are imposed on each iteration, where we take Rh = 30 kpc. 

For a given Zh the diffusion coefficient as a function of momentum and the reacceleration or convection 
parameters is determined by boron-to-carbon (B/C) ratio data. The spatial diffusion coefficient is taken 
as D xx = (3Dq{p I po) s if necessary with a break (5 = 5± below rigidity po, 5 = 82 above rigidity po). The 
injection spectrum of nucleons is assumed to be a power law in momentum, dq(p)/dp oc for the injected 
particle density. For the ca se of reacceleration the momentum-space diffus ion coefficient D pp is related to 



the spatial coefficient D xx ( Berezinskii et all . Il99dl : Seo and Ptuskinl 19941 ) via the Alfven speed va- The 
convection velocity (in z-direction only) V(z) is assumed to increase linearly with distance from the plane 
(dV/dz > for all z). This implies a constant adiabatic energy loss. 

The interstellar hydrogen distribution uses H 1 and CO surveys and information on the ionized component; 
the helium fraction of the gas is taken as 0.11 by number. The H2 and H 1 gas number densities in the 
Galactic plane are defined in the form of tables, which are interpolated linearly. The extension of the gas 
distribution to an arbitrary height above the plane is made using analytical approximations. 

The distribution of CR sources is chosen t o reproduce the CR d istribution determined by ana lysis of 
EGRET 7-ray data ( Strong and Mattoxl . 199fil ) and was described in Strong and Moskalenkol (|l998l ). 

Energy losses for nucleons by ionization and Coulomb interactions are included, and for electrons by 
ionization, Coulomb interactions, bremsstrahlung, inverse Compton, and synchrotron. 

Positrons and electrons (including secondary electrons) are propa gated in the same model. Positron 
production is computed as described in Moskalenko and Strongl |l998), that paper includes a critical reeval- 
uation of the secondary tt^- and .fT^-meson decay calculations. 
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Gas-related 7-ray intensities are computed from the emissivities as a function of (R,z,E^) using the 
column densities of H 1 and H2. The interstellar radiation field, used for calculation of the inverse Compton 
emission and electron energy losses, is calculated based on stellar population models and COBE results, 
plus the cosmic microwave background. 

New Developments 

The experience gained from the original fortran-90 code allowed us to design a new version of the model, 
entirely rewritten in C++, that is much more flexible. It allows essential optimizations in comparison to 
the older model and a full 3-dimensional spatial grid. It is now possible to explicitly solve the full nuclear 
reaction network on a spatially resolved grid. The code can thus serve as a complete substitute for the 
conventional "leaky-box" or "weighted-slab" propagation models usually employed, giving many advantages 
such as the correct treatment of radioactive nuclei, realistic gas and source distributions etc. It also allows 
stochastic SNR sources to be included. It still contains an option to switch to the fast running cylindrically 
symmetrical model which is sufficient for many applications such as the present one. 

In the new version, we have updated the cross-section code to include the latest measurements and energy 
dependent fitting functions. The nuclear reaction network is built using the Nuclear Data Sheets. The 
isotopic cross section database consists of thousands of points collected from the literature. This includes 
a critical re-evaluation of some data and cross checks. The isotopic cross sections are calculated using the 
author's fits to major beryllium and boron production cross s ections p + C, N, O — » Be, B. Other cross sec- 
tions a re calculated using phenomenological approximations by Webber et al. ( 1990h and/or lsilberberg et al 



1998) renormalize d to the data where it exists. The cross sections on the He target are calculated using a 



parametrization bv lFerrando et al.l 

The reaction network is solved starting at the heaviest nuclei (i.e., 64 Ni). The propagation equation 
Eq. (^Q) is solved, computing all the resulting secondary source functions, and then proceeds to the nuclei 
with A — 1. The procedure is repeated down to A = 1. In this way all secondary, tertiary etc. reactions 
are automatically accounted for. To be completely accurate for all isotopes, e.g., for some rare cases of 
/3^-decay, the whole loop is re peated twice. Our preliminary results for all CR species Z < 28 are given in 



Strong and Moskalenknl fcOQlh. 



COSMIC RAYS IN INTERSTELLAR TURBULENCE: A SELF-CONSISTENT APPROACH 

It is known that galactic CR have relatively high energy d ensity and they can not always be treated as 



test particles moving in given interstellar magnetic fields, see Be rezinskii et a l. (1990). In particular, when 
the stochastic reacceleration of CR is considered, one has to take into account the damping of waves, which 
loose energy on particle acceleration. The damping causes the change of the wave spectrum that in its turn 
affects the particle transport. Thus in principle the study of CR propagation may need a self-consistent 
consideration. We shall see from the numerical estimates below that the back reaction of CR on interstellar 
turbulence should be taken into account at particle energies below 10 GeV/nucleon. 

It is assumed that the wave-particle interaction in the interstellar plasma has a resonant character (the 
cyclotron resonance). The scattering of energetic particles by random hydromagnetic waves leads to the 
spatial dif fusion of CR in the Gala xy. The CR diffusion coefficient is determined by the following approximate 
equation ( Berezinskii et al. . 1990h : 

D{p) = i27rk re 9 s w(k res y (2) 

where v is the particle velocity, r g = pcj (ZeB) is the particle Larmor radius in the average magnetic field 
B (Ze is the particle charge), k res = l/r g is the resonant wave number, and W(k) is the spectral energy 
density of waves defined as J dkW(k) = 5B 2 /4ir (5B is the random magnetic field, 5B <C B). 

In its simplified form, the steady state equation for hydromagnetic waves with a nonlinear energy transfer 
in /c-space can be written as 

^(c A k 2 W(k)J^^\ = -2T cr (k)W(k) + SS(k-k L ), (3) 
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k > ki (e.g., Landau and Lifschit d . 19871 : Norman and Ferraral . 1996). Here the l.h.s. of equation describes 
the Kolmogorov type nonlinear cascade from small k to large fc, C a is a constant and according to the 
numerical simulations of magnetohydrodynamical turbulence by Verma et al. ( 1996h the magnitude of Ca 
is very roughly equal to 0.3, p is the interstellar gas density. The r.h.s. of Eq. (j3J) includes the wave damping 
on CR and the source term which works on a large scale ~ 1/&l and describes the generation of turbulence 
by supernova bursts, powerful stellar winds, and superbubbles expansion. In the limit of negligible damping, 
r cr = 0, the solution of Eq. (j3J) gives the Kolmogorov spectrum W(k) oc /c~ 5 / 3 at k > kj_,. The latter leads 
to the diffusion coefficient D(p) oc v(p/Z) 1 ^ 3 . 



The expression for wave amplitude attenuation is given by (|Berezinskii et al. , 1990): 



r, 



where V a = B /yJAnp is the Alfven velocity, p res 
The solution of Eqs. (JHJ), (jlj) at k > kz gives 
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ZeB/ck is the resonant momentum. 
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The wave damping on CR is essential only at relatively small scales k 1 < k d 1 , k^ <C k L 1 (in fact, k d 1 ~ 10 
cm, k L x ~ 10 20 cm). The turbulence at k d 1 <C k 1 < k L x obeys the Kolmogorov scaling W{k) oc k 5 / 3 . It 
is clear from Eq. © that the dissipation of waves on CR described by the second term in square brackets 
makes the wave spectrum at large wave numbers steeper than the Kolmogorov spectrum. 
The diffusion mean free path l(p), as defined by D(p) = vl(p)/3, is now given by: 
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where we use the resonant conditions p = p re s(k) and introduce pi = Presiki)- Eq. © is formally valid 
at p < pl ~ 10 17 eV/c. The second term in square brackets reflects the modification of diffusion mean 
free path due to the wave dissipation on CR. It increases with decreasing the particle energy. It is easy to 
estimate that the second term in square brackets is small at high energies but it is comparable with the 
first term in square brackets at p ~ 1 GeV/c (at / = 3 pc, V a = 30 km s _1 , B = 3 pG). So, the nonlinear 
modification of CR diffusion is expected at E < 10 GeV/nucleon. 

As the most abundant specie, the CR protons mainly determine the wave dissipation. Their distribution 
function ^(p) should be used to calculate W(k) and l(p). The self-consistent treatment of CR propagation 
implies the solution of transport Eq. Q for ^(p) with the diffusion mean free path Eq. ©, which is a 
function of ^(p). 

To demonstrate the effect of wave damping, let us consider a simple case of one-dimensional diffusion 
model of CR propagation with the source distribution q(r,p) = qo(p)5(z) (that cor responds to the infinitel y 
thin disk of CR so urces located at z 
( 200lh . 



Jones et al 



0) and the flat CR halo of height H, see iBerezinskii et al. 
The CR source spectrum is of the form qo(p) oc p~~ /s and the index is estimated as 
7 S = 2.0. ..2.4. 

Let us assume that stochastic reacceleration does not essentiality change the CR spectrum on the charac- 
teristic time of CR leakage from the Galaxy. We also ignore ionization energy losses and losses for nuclear 
fragmantation and assume that the equilibrium density of CR protons is determined by the balance between 
their production in the sources and the free diffusion leakage from the Galaxy. The solution of the diffusion 
equation for the relativistic protons in the galactic disk is then 



(7) 
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We simplify Eq. © by using the approximation ^■^(pi) ~ ^{tpi) and write it down as 



l(p)~l K (p) 



(8) 



where Ik(p) 1S the diffusion mean free path in the turbulence with Kolmogorov spectrum, which is not 
modified by the wave damping on CR. 

Now one can obtain the following solution of Eqs. (JJJ) and © for the diffusion mean free path in the case 
of a power law source spectrum: 



The last expression together with Eq. Q presents the self-consistent solution of the considered simple 
problem. 

The low-energy asymptotic mean free path Eq. Q is l(p) oc v~ 2 p 3 ~~ 2,ys . The high-energy asymptotic mean 
free path Eq. is l(p) ~ Ik(p) P 1/ ' 3 that correspons to a Kolmogorov spectrum of waves. One can check 
that the wave spectrum preserves the Kolmogorov scaling W(k) oc /c -5 / 3 at small k and it is modified as 
W(k) oc k l ~ 2ls at large wave numbers. 

CONCLUSION 

The back reaction of energetic particles on interstellar turbulence may change the spectrum of hydromag- 
netic waves at large wave numbers that leads to the change of CR diffusion coefficient. In particular, the 
modification of Kolmogorov type spectrum of turbulence leads to the CR diffusion coefficient which goes 
through a minimum at particle magnetic rigidity about a few GV, so that the diffusion coefficient has the 
power law asymptotics D(p) oc vp 1 / 3 and D(p) oc v^ 1 p 3 ~ 2 ^ /3 at large and small p respectively (the exponent 
of CR source spectrum 7 S ~ 2.2). Such a behaviour reduces the rate of stochastic reacceleration of CR in the 
interstellar medium that is needed to reproduce the observerd peaks in the secondary-to-primary elements 
ratios in CR at these rigidities. The reduced rate of reacceleration may help to explain the antiproton data. 

In a future work, we plan to undertake a full scale self-consistent modeling of CR transport in the turbulent 
interstellar medium in the frameworks of the realistic galactic model described here and with the damping 
of waves on energetic particles included. 
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